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We study the origin of the structure in the spin-polarized 
electron energy loss spectroscopy (SPEELS) spectra of ferro- 
magnetic crystals. Our study is based on a 3d tight-binding 
Fe model, with constant onsite Coulomb repulsion U between 
electrons of opposite spin. We find it is not the total density 
O , of Stoner states as a function of energy loss which determines 
the response of the system in the Stoner region, as usually 
thought, but the densities of Stoner states for only a few in- 
. terband transitions. Which transitions are important depends 
ultimately on how strongly umklapp processes couple the cor- 
! t | responding bands. This allows us to show, in particular, that 
■ the Stoner peak in SPEELS spectra does not necessarily in- 
dicate the value of the exchange splitting energy. Thus, the 
I , common assumption that this peak allows us to estimate the 
magnetic moment through its correlation with exchange split- 
ting should be reconsidered, both in bulk and surface studies. 
Furthermore, we are able to show that the above mechanism 
is one of the main causes for the typical broadness of experi- 
mental spectra. Finally, our model predicts that optical spin 
waves should be excited in SPEELS experiments. 
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I. INTRODUCTION 

The study of elementary excitations in itinerant- 
electron ferromagnets is an area which is currently very 
active in spite of the enormous amount of publications on 
the subject since early work in the sixties (see e.g. Refs. 
jl],^) and references therein). From the beginning, exper- 
imental and theoretical work on these materials concen- 
trated on neutron scattering and dynamical susceptibil- 
ity studies. (l]|||6| Efforts have continued in this direc- 
tion until today, both because of the gradual improve- 
ment of electronic band structure calculations [0J|] and 
because of the improvement of the experimental method. 
H Around the mid-eighties, however, a new technique 
was introduced in the field, namely, spin polarized elec- 
tron energy loss spectroscopy (SPEELS). Among its first 
successes, one can count the first observations, in a fer- 
romagnetic glass and in nickel, (ll]] of what were in- 
terpreted as Stoner excitations. Further work, reporting 
more detailed measurements, confirmed those findings. 
p2| , p^| Theoretical model calculations of inelastic elec- 
tron spin-flip exchange scattering, Jl^-|r^| provided a ba- 
sis for the interpretation of those observations in terms 
of Stoner excitations. In addition, Vignale and Singwi 
found in their work that spin waves should also be ob- 
servable in SPEELS measurements. jl6| However, these 
had not been observed at the time, nor were they ob- 
served in the several years that followed. Spin waves 
were found in other model calculations, (l^Jl^l the cal- 
culations of Plihal and Mills being the most conclusive 
in this respect because of their more accurate treatment 
of electronic structure. jn| It is only very recently that 
the detection of spin waves in a SPEELS experiment has 
finally been reported. Q The application of SPEELS 
has been naturally extended to the study of magnetic 
surfaces. pij-p3[ An important theoretical effort in this 
direction is that by Mills and collaborators, p^-[2^| who 
have studied ferromagnetic thin films as well. 

To introduce the questions addressed by this work, 
we recall briefly some of the main concepts involved 
in SPEELS and discuss some of the findings to date. 
SPEELS is a spin-polarized version of electron energy 
loss spectroscopy in the sense that the spin polarization 
of the scattered electrons is also measured. The imping- 
ing electron often is also spin-polarized, but this is not 
necessarily so (see e.g. Refs. In a so-called spin- 

flip exchange scattering event, an incoming electron with 
given spin comes to occupy an empty level in the material 
while an electron with opposite spin is driven out and is 
detected. The process produces thus a Stoner excitation. 
In the band picture of magnetic transition metals, be- 
low the Curie temperature, the exchange split 3c? bands 
provide large densities of occupied majority-spin states 
below the Fermi energy and vacant minority-spin states 
above it. Thus, it is more likely for an impinging elec- 
tron with minority spin to excite a Stoner pair than for a 
majority-spin incoming electron, particularly for an ex- 



citation energy corresponding to the exchange splitting 
of the ferromagnct. This is the mechanism invoked to 
explain the Stoner peak or the asymmetry reported in 
Ref s. [[Tof and and further experimental work (Refs. 
Jl2| , |l3[ ). However, it turned out necessary to elaborate 
on several other issues. Firstly, the Stoner peak was 
very broad in all observations. This was interpreted by 
Kirschner, Rebenstorff, and Ibach as an indication of 
the nonuniformity of exchange splitting throughout the 
Brillouin zone. Then, in Fe, the energy loss at which the 
Stoner peak occurs and its width were reported by Venus 
and Kirschner to increase with increasing scattering an- 
gle, jl2| a fact that was correlated by these authors with 
the calculated density of Stoner states. Also, a thresh- 
old for the onset of Stoner excitations in Ni(llO) was 
reported by Abraham and Hopster | fl3| and interpreted 
in terms of the Ni 3d band structure. These workers, 
moreover, indicated that their spectra did not differ sig- 
nificantly for off specular scattering angles ranging from 
10° to 40°, [^7| which they explained as due to the non- 
conservation of the momentum component perpendicular 
to the surface. 

Finally, an important application based on SPEELS in- 
terpretation is that, in surface and thin film studies, the 
Stoner peak is assumed to give information on the sur- 
face magnetic moment through the correlation between 
exchange splitting and moment. |^,^3| In particular, a 
Stoner peak found at higher energies than the exchange 
splitting bulk value is assumed to indicate an enhanced 
magnetic moment at the surface. 

Clearly, more theoretical work is required, for the bulk 
as much as for surfaces, to make further progress. In par- 
ticular, it would be important to understand better the 
phenomenology of SPEELS and to try to be more specific 
about the information we can expect from it. This would 
also provide experimenters with useful feedback. Accord- 
ingly, we think it is worthwhile going back to a model cal- 
culation and look more closely at the dynamic properties 
of the material probed by SPEELS. In this work we con- 
sider a model of Fe based on paramagnetic tight-binding 
3d bands, with up and down spin bands rigidly split. The 
cross section for spin-flip exchange scattering processes 
is evaluated within the random phase approximation, as- 
suming the solid is described by a multi-band Hubbard 
Hamiltonian. As we shall see, this allows us to show that 
it is not the total density of Stoner states as a function 
of energy loss and momentum transfer which causes the 
structure in the Stoner region of the spectrum, but the 
density of Stoner states for a few interband excitations. 
Which interband excitations are important is essentially 
determined by the weight of the matrix elements for such 
processes. It this regard the contribution of umklapp 
scattering is fundamental because of the the coupling of 
different bands at different energy ranges. This gives rise 
to a richer structure in the Stoner region of the spectra. 
Also, our model predicts that optical spin waves should 
be observable through SPEELS. Again umklapp scatter- 
ing proves critical, providing optical spin waves with the 
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necessary oscillator strength. We find this strength to 
depend importantly on scattering angle. 

Section II of this paper is devoted to theory, presenting 
the derivation of the spin-flip exchange scattering cross 
section for our model. We present our main results in 
Section III. Then follows in Section IV a discussion of 
our results in the light of experimental findings and other 
theoretical work. Finally, in Section V we summarize our 
work and give some conclusions. 

II. THEORY 

A. Spin-flip exchange scattering cross section 

The electron spin-flip exchange scattering differential 
cross section for a iV-electron system target has been 
previously derived on general grounds by Vignale and 
Singwi in terms of a particle-hole excitation correlation 
function. jl6| One has 

(Pa m 2 p f 1 Imx£U f (P»,q,ff) 

dEdQ~ 47r 2 ?i 4 p 4 7r l-e-P E ' ^ ' ' 

where E is the energy loss, f2 is the solid angle, m the 
electron mass, and /3 = l/fc^T. Momentum transfer is 
given by q = pi — p/ , with pt and p / the momentum of 
the incoming and outgoing electrons, respectively. Like- 
wise, ai is the spin of the impinging electron and cr/ that 
of the scattered one. The retarded function X^o-* can 
be obtained by analytic continuation of the two-particle 
temperature correlation function 

f fi - 

Xo-*07(pi,q, ioJ n )= ~ / dre %UJnT 
Jo 

x ( T T [g aia} (jpi, q, T) e l. af (pi, q)] ). (2.2) 

In this equation lu u — 2irn/ (3 is a bosonic Matsubara fre- 
quency, T T is the imaginary time ordering operator, and 
the brackets indicate the thermodynamic average in the 
canonical ensemble. [|8| £>£. i(T/ is the particle-hole creation 
operator 

N 

eU,(p*>q) = "jvE / llv,lv <' r <'" r 

3=1 J 

xvdr-T^Kr)^^), (2.3) 

where V4-( r ) is the field operator creating an electron of 
spin a at position r and v(r) — e 2 /r is the Coulomb 
interaction between the scattered and target electrons. 
The sum runs over the N electrons in the target system. 
This expression is quite general, and could be applied 
equally well to a solid, an atom, or a molecule. 

We now consider the N electrons in a crystal mate- 
rial. We write the Bloch wave function for a state with 
wave vector k and spin a in band n in terms of Wannier 
functions, 



^ nk(T (r) = 4rE e lk ' R 0„ k (r - R)r, CT , (2.4) 
ViVo R 

where A*o is the number of sites in the crystal and r\ a is 
the spin function. Denoting by ajj kff the operator creat- 
ing an electron in such a state, the field operators can be 
expanded as ipl(r) = J2 n k Vw(r)o4 k(T . The particle- 
hole creation operator becomes 

rin' k 

(2.5) 

where the sum in momentum space runs over the Bril- 
louin zone and matrix element W nn i is given by 

WW(p,,q,k) = ^^f)(k- Pi -K)0* k (k-K) 

K 

x^ k - q (k-q-K). (2.6) 

Here K denotes vectors in the reciprocal lattice and V 
is the volume of the sample. The * indicates a Fourier 
transformed function and * denotes complex conjuga- 
tion. To write the last two equations we have defined 
a nk+Ka = a nker ana - nave exploited the periodicity of 
the Wannier functions in the wave vector index, i.e. 

0nk+K<7 = </>nk<T- 

B. RPA expression for a tight-binding system 

We are interested in the cross section for an itinerant 
electron ferromagnet. We describe the system within a 
tight-binding approximation, thus writing the Wannier 
wave function for given k and band index n by a linear 
combination of atomic orbitals ip m 

^nk(r) = y^bmn(k)y w (r)- (2.7) 

m 

The coefficients b mn diagonalize the crystal Hamiltonian 
and are normalized so as to define a unitary matrix. 
Consequently, the independent-electron Hamiltonian of 
the system can be written Hq — e n (p)aJ l p (T (i n p (T , 
where the e n (p) are paramagnetic band energies. We as- 
sume the interacting system is described by a multi-band 
Hubbard Hamiltonian. In our basis we have 

fff = ^ I! I! c «™(p + q',p)<vm'(p' - q',p') 
° nn '' a pp/q/ 

Xa np+q'ff a »»P<7 a n'p'-q'-(T a m'p'-o'> (2-8) 

where we have defined c„ m (p, q) = J2i &in(p)&Zm(q), and 
U is the effective on-site Coulomb interaction for two 
electrons with opposite spin. The RPA evaluation of the 
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correlation function Xo- i0 -. defined in Eq. 2.2 is a straight- 
forward generalization of that in previous work fisfl . The 
response function divides naturally in two, 



(2.9) 



The Stoner or single-particle excitation contribution is 
given by 



x: 



/n'k— 



qcr f 



f, 



nn' k 



e n / CT/ (k - q) - e nrTi (k) 



x|W nn /(pi,q,k)| : 



(2.10) 



We have introduced the occupation probability of state 
nka, f n ka = (a* kcr a„ k(T ), and the single-particle energy 
modified by the exchange self-energy 



e nCT (k) = e„(k) 



U 
An 



f npcr ■ 



(2.11) 



Thus, in this model, spin down and spin up energy bands 
are rigidly split by the quantity A = U({n^) — (n^)), 
where 



N 



(2.12) 



is the average number per site of states with spin a. 
The many-body contribution is given by 



U 



nn' 

xr™'( Pi ,q,iu, n ), (2.13) 



with the auxiliary functions G nn and T nn defined as 
follows: 

s ^ l / — d iuj n + e m , a (k - q) - e mai (k) 

mm' k 1 

xWC m ,(pi,q,k)& nm (k)b n / m ,(k- q), (2.14) 

and, considering G and Y as vectors with coefficients in- 
dexed by nn', 



U 



1 + TrB<r (ff/ (q,iWn) 



r<7 jCT/ (Pi,q,iw n ) 

x G CTlCr/ (p i ,q, iuj n ] 
where the elements of matrix D are 
D^ f mm '(q,iu; n ) = ^E^(k)^(k- q) 



(2.15) 



W k 

/i'k— qcr/ fi\i.o~i 
iu n + tva; (k - q) - e; CTi (k) 



fe m; (k)6 m - r (k-q). (2.16) 



To obtain X^ ia p analytic continuation ioj n — > E + ir\ of 
the above results is straightforward. Also, results in next 
section are obtained in the zcro-tcmpcraturc limit. 



III. RESULTS 

Before presenting our results, there are a few details of 
our model we should discuss. As we said in the introduc- 
tion, our Fe model is based on simple 3d tight-binding 
paramagnetic bands, i.e. we neglect sp hybridization. 
This has, of course, incidence on the quantitative de- 
tails of the response of the system. However, we will 
see our results are rather general and do not depend on 
these details. The Wannier functions were written as a 
linear combination of the Fe five 3d atomic wave func- 
tions. The overlap integrals were calculated, up to next 
nearest neighbors, using the Fe atomic wave functions 
determined according to Griffith's prescription, p9[ and 
a lattice constant a — 2.87 A. We have, thus, a five 
band model. The bandwidth was set to 4.7 eV, which 
corresponds roughly to the bandwidth of d electrons in 
Fe. 1 31 1 Exchange splitting A was chosen to be 2 eV, 
taking as reference the position of the peaks in the den- 
sities of states for up and down spins. Then, the Fermi 
level was fixed by the condition of having six electrons 
per unit cell. We show the density of states for up and 
down spins in Fig. |l|. These exhibit the bonding and anti- 
bonding regions common to BCC materials with unfilled 
d shells. Once the Fermi energy is fixed we can deduce 
the strength of the effective Coulomb interaction U in 
our model from U — A/((n|) — (^x))- We find U = 0.69 
eV, which compares very well with the energy found in 
other works. |32| On the other hand, bulk polarization 
is too high, roughly 48%, (33| reflecting the lack of hy- 
bridization with sp electrons. 
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FIG. 1. The Fe up and down spin electron densities of 
states for our model, exhibiting the characteristic bonding 
and antibonding regions. Bandwidth is 4.7 eV and exchange 
splitting 2 eV. 
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In a typical SPEELS experiment, the incoming elec- 
tron beam impinges on the sample surface at an angle 
9 to the normal and the total scattering angle is 90°. 
O |l3| , p2|j23| For fixed scattering angle, i.e. for given in- 
coming and outgoing momenta, there are three possible 
scattering processes corresponding to different momen- 
tum transfer. In two relatively small angle in- 
elastic scattering event is preceded or followed by elastic 
scattering. In the third one, all the momentum transfer 
is absorbed by the electron-hole pair excitations, i.e. it 
is a large angle scattering event. We consider here the 
geometry of Venus and Kirschner, jl2) that is, the sam- 
ple exposes the (110) surface and the scattering plane is 
defined by the surface normal [110] and the [001] axis. If 
u is the axis normal to the surface, the three momenta 
mentioned are given as a function of energy loss E and 
impinging momentum pi and energy Ei by 



Total ■ 

q. q' 



q u = Pi (cos 9 - sin 6y/l - E/Ei), 
q' u = Pi (-cos0 + &m0y/l-E/Ei), 



-Pi (cos 6 + sin 6y/l - E/Ei), 



Pi(sin9 - cos 6y/l - E/Ei), 



(3.1) 



where 9 is the angle to the normal. The large angle scat- 
tering event corresponds to q". Momentum transfer par- 
allel to the surface is the same in the three cases. To cal- 
culate the final spectrum, the contributions of these three 
processes have to be added because experiment does not 
discriminate them. 

Also, since Fe presents a non negligible quantity of free- 
like s and p states at the Fermi surface, the interaction 
between the incoming electrons and those in the solid 
will be screened. We take this into account using the 
Thomas-Fermi form of the screened Coulomb potential, 
with a screening wave vector corresponding to the density 
of states of s and p electrons at the Fermi surface. |34|] 
This gives (/tf = 0.26 in units of k a = An /a. Finally, an 
important point to mention is that we use a finite value 
for rj when taking the analytic continuation ioj n — > E+irj. 
Since, to our knowledge, there are no estimates of the self- 
energy corrections for up and down spin bands in Fe, ]35f] 
we take rj — 80 meV, which corresponds to the resolution 
in the latest experiment on this material. pOfl 



A. Interband densities of Stoner states 

Let us consider a majority spin electron with an angle 
of incidence 9 — 60° to the normal and an incoming en- 
ergy Ei = 22 eV, which is the energy used in Ref. Jl2| . 
We see in Fig. ^(a) that the total spin-flip exchange scat- 
tering cross section is indeed rather broad, with its peak 
centered at an energy much higher than the exchange 
splitting value (2 eV in our model), a trend observed ex- 
perimentally by Venus and Kirschner [Q . We also show 
the partial cross sections for different momentum trans- 
fers. The curves for small scattering angle coincide for 
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FIG. 2. (a) Spin-flip exchange cross section for a majority 
spin electron, with impinging energy of 22 eV and angle of 
incidence of 60°. We show the total cross section, as well 
as the partial cross sections for different momentum transfer. 
The peak is at 3 eV, an energy much higher than the exchange 
splitting A = 2 eV. (b) Interacting and noninteracting cross 
sections. The difference between both curves clearly shows 
two collective modes, one just below 2 eV, and the other below 
1 eV. The noninteracting cross section shows three distinct 
features, namely, the peak at 3 eV, a shoulder at 4 eV, and a 
broad hump around 1 eV. (c) The total and partial densities 
of Stoner states. These show the typical maxima at exchange 
splitting, which are absent from the SPEELS spectrum. The 
densities of Stoner states are incapable of explaining the peak 
of the SPEELS spectrum at 3 eV. 

symmetry reasons. Though total cross section broadness 
is somewhat increased because of the difference between 
small angle and large angle scattering, the cross section 
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in each case is broad in itself. We have examined the 
origin of the structure in this spectrum. Firstly, we sep- 
arated single-particle excitations and many-body effects. 
In Fig. ^(b) we show the noninteracting and interacting 
(total) cross sections. This figure clearly shows us the 
contributions of collective modes. Indeed, the broad fea- 
ture starting around 0.2 eV indicates the excitation of low 
lying spin waves, and, more interestingly, the shoulder at 
higher energy, below 2 eV, indicates the excitation of op- 
tical spin waves. This is important because optical spin 
waves have not been discussed previously in connection 
with SPEELS measurements. We consider spin waves 
again further on and we concentrate here on the single 
particle traits. Besides the peak at 3 eV, the noninter- 
acting cross section shows a shoulder at 4 eV and a broad 
feature, albeit much smaller, at low energy loss, around 
1 eV or so. SPEELS spectra have often been interpreted 
in terms of the density of Stoner states. Accordingly, we 
show in Fig. ||(c) the total density of Stoner states, as well 
as the densities of Stoner states for different momentum 
transfer (as before, the curves for small angle scattering 
are the same) . |3(| It is evident in the cross sections that 
there is nothing reminiscent of the high density of states 
at the exchange splitting energy. The only features of the 
cross sections that can find an explanation in the density 
of Stoner states are the shoulder at 4 eV and, possibly, 
the hump around 1 eV. We have, thus, refined our study 
and have considered the behavior of the density of Stoner 
states as a function of energy loss and the bands coupled 
in an excitation (recall energy loss and momentum trans- 
fer are coupled, cf. Eq. (3.1)) 



^nn'iE) — ^ ] (/n'k— qg-f /rikffi ) 
k 

xS{E + e nV , (k - q) - e„ CTj (k)). (3.2) 



Hence, subscripts n and nl indicate minority and major- 
ity bands, respectively (bands are numbered from bot- 
tom to top). We show a plot of the density of states thus 
defined in Fig. ||(a). We can see that the allowed and 
forbidden interband excitations are completely identified 
(we use the term forbidden for interband excitations with 
vanishing density of Stoner states). Moreover, the series 
of Stoner peaks clearly reflect the bonding and antibond- 
ing nature of the electronic structure, giving rise to two 
arrays of peaks, for higher and lower excitation energies. 
p7| The question is, of course, which of these Stoner 
peaks contribute the most to SPEELS cross sections. 
The answer is to be found, perhaps unsurprisingly, in 
how strongly the different bands are coupled by the ma- 
trix elements W nn i of the electron-hole creation operator, 
which weights the contribution of each Stoner excitation 
(see Eq. 2.10| ). What is not so obvious is the outcome of 
the combined effect of Stoner peaks and matrix elements. 
Let us consider the average of the square of the absolute 
value of matrix elements W nn > over the Brillouin zone, 
i.e. 
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FIG. 3. (a) 3D plot of the interband densities of Stoner 
states. Allowed and forbidden interband excitations are 
clearly identified. The two series of peaks reflect the bonding 
and antibonding nature of the Fe electronic structure. The 
highest peak is found at exchange splitting, for nn' — 44. 
As we explain in the text, however, this peak contributes lit- 
tle to the spin-flip exchange cross section, (b) 3D plot of 
the average ( | W nn < j 2 ). The slight dependence on energy re- 
sults in the wave like form of the surface. The lines on the 
surface parallel to the energy axis correspond to fixed band 
couple value. The important interband excitations are deter- 
mined by the crests, namely 21,22,31,32,41, and 42. One can 
also observe that ( |Wn n '| 2 ) reaches its lowest values for n5, 
with n — 1, ... ,5. (c) 3D plot of ( \ W nm \ 2 ) without umklapp 
processes, The important interband excitations have been re- 
duced to nn' = 31, 32, thus singling out excitations in a re- 
stricted energy range. 



\w n 



(3.3) 



G 



(v denoting the volume of the Brillouin zone). In 
Fig. ||(b) we show the graph of ( |W„„'| 2 ) as a function 
of energy loss and of bands coupled. There is little sig- 
nificant variation as a function of energy loss, but a very 
important structure as a function of band couple, result- 
ing in a wave like pattern. Comparing Figs. §(a) and [§(b) 
we can clearly see when it is that both quantities, p nn > 
and (|W rm '| 2 ), interfere constructively. Thus, although 
the density of Stoner states reaches is highest peak at ex- 
change splitting, the average ( IW-WI 2 ) is negligible for 
the corresponding band couples. Instead, although the 
densities of Stoner states for interband excitations 21 and 
22, 31 and 32, and 41 and 42 are more modest, the cor- 
responding matrix element averages are high. Whence 
the peak around 3 eV and the shoulder around 4 eV in 
the noninteracting cross section in Fig. ^(b). Actually, 
the peak at 3 eV is more of a hat on top of the high 
cross section value due to interband excitations 31 and 
32. We can also see that the hump around 1 eV is due 
to excitations coupling bands 2 and 3, and 2 and 4. To 
corroborate our analysis, we show in Fig. || the cross sec- 
tion taking into account solely the interband processes 
mentioned above. We include the total noninteracting 
cross section for comparison as well. We see that the 
few interband excitations considered indeed account al- 
most completely for the structure of the noninteracting 
spectrum. An argument to understand how so simple a 
picture can work is that, since the atomic 3d orbitals are 
localized, their Fourier transform is rather flat, so that 
\W nn '\ 2 in the single-particle correlation function Xa at 
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No umklapp: int 
No umklapp: nonint. 



(cf. Eq. ( 2.10 )) may be replaced by its average value over 
the Brillouin zone. Thus, X%af * s approximately propor- 
tional to J2nn' Pnn'(\W nn >\ 2 ), a weighted average of the 
interband densities of Stoner states. 
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FIG. 4. Noninteracting scattering cross section taking into 
account only 8 (W = 21,22,23,24,31,32,41,42) out of the 
25 possible interband excitations, selected as explained in the 
text, compared to the total noninteracting scattering cross 
section. The first curve follows very closely to the second 
one. 
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FIG. 5. Spin-flip exchange scattering cross section with- 
out umklapp processes. All that is left is a broad maximum 
around 3 and 4 eV, and the small hump at small energy loss. 
Clearly, the most important information is lost, i.e. the peak 
at 3 eV and the shoulder at 4 eV (cf. Fig. |(b).) 



B. Umklapp processes 

Another most interesting phenomenon playing a fun- 
damental role in SPEELS is umklapp scattering. One 
can see in Eq. (2^) that the contribution of umklapp 



processes to the particle- hole excitation operator g ai<Jf is 
weighted by the Coulomb interaction and the Wannier 
wave functions. Because of the decay of the Coulomb 
potential as well as of the atomic orbitals with increas- 
ing wave vector, the weight becomes rapidly negligible 
for lattice vectors beyond first nearest neighbors. This is 
enough, however, for umklapp processes to have a twofold 
effect. To see this, let us consider cross sections taking 
into account only normal excitations (i.e. with respect to 
the first Brillouin zone). We show this in Fig. [|, where 
we plot both the interacting and noninteracting cross sec- 
tions. Firstly, we see that, quite apart from their much 
lower values in comparison with the full response case, 
spectra in Fig. ^| show little resemblance with those in 
Fig. ||(b) (scale in both figures is the same). This is be- 
cause the possible interband excitations have been drasti- 
cally reduced. Indeed, let us consider the graph of the av- 
erage ( | W nn i \ 2 ) for excitations strictly conserving crystal 
momentum. We show this in Fig. ||(c). The wave crests 
have been reduced to that for nn' = 31,32, from which 
it is obvious the the different wave crests in Fig. |](b) are 
due to umklapp scattering. Normal scattering alone re- 
sults in a spectrum almost completely distorted because 
of the excitation of interband transitions mainly for en- 
ergies between 3 and 4 eV (cf. Fig. ||(a)). 

Also, in Fig. || it is immediately apparent that there 
remains no trace of spin waves in the spectrum. Indeed, 
the interacting and noninteracting curves are almost in- 
distinguishable, with no hint of the spin wave modes be- 
low 2 or 1 eV. Thus, it appears that umklapp processes 
provide collective excitations with oscillator strength, at 
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FIG. 6. Contour plot of (IM 7 ^'! 2 ) as a function of scat- 
tering angle, for an energy loss E = 2.5 eV. Thinner line 
correspond to higher averages. Drastic changes occur partic- 
ularly for near specular scattering. In this case, the interband 
excitations for nn' — 22, 33 and 44 will overwhelmingly dom- 
inate the spectrum. Averages decrease significantly for large 
scattering angles. Comparing values for 60° and for 55°, we 
see that interband excitations for nn' — 41,42 are less im- 
portant for the latter than for the former. Similarly, for 70°, 
interband excitations 21 and 22 play a less important role 
than for 60°. 

least optical modes. We must point out here that the 
low lying collective mode contributing to the broad fea- 
ture below 1 eV is not an acoustic mode, but also an 
optical one. We have calculated its dispersion relation 
in the [100] direction and found that it tends linearly to 
218.7 meV for q — ► 0. The slope is positive, but very 
low, with a value of 8.9 meV- A. The reason is that the 
energy bands in our model are purely d. It is well known 
that models of itinerant ferromagnetism which do not 
take into account hydridization with sp bands, fail to 
describe spin waves appropriately. |Q So we do not ex- 
pect our model to predict accurate dispersion relations 
for collective excitations. However, we do think our re- 
sult properly introduces optical spin waves as a source of 
structure in SPEELS measurements. 



C. Spectra and angle of incidence 

To illustrate further the pertinence of our analysis, let 
us consider an example. A question addressed in the past 
by experimenters, without finding a clear answer, [ ^2|Ji~3| 
has been that of the variation of spectra with scattering 
angle. Since (|W„„<| 2 ) plays such a consequential role, 
we take a look at its dependence on angle of incidence 
through the contour plot in Fig. ||. Generally speak- 
ing, the same interband averages remain the most im- 
portant important as angle changes, except toward spec- 
ular scattering, when weights raise and shift significantly 
(then the highest values are to be found for nn' = 22, 33, 
and 44). Let us consider the spectra for angles of inci- 
dence 9 = 55° and 9 = 70°. In the first case, Fig. | 
shows us that the importance of interband excitations 



nn' = 41,42 is diminished with respect to the 9 = 60° 
case. For 9 — 70°, it is band couples nn' = 21, 22 that are 
diminished. In this way, in fact, we hinder the contribu- 
tion of umklapp processes coupling different bands. Re- 
garding the densities of Stoner states, on the other hand, 
we found that changes from one angle to another are 
rather small, affecting essentially only the height of the 
peaks. This is because a change in angle will not change 
the energies at which there can be a Stoner excitation, 
but basically the number of these. Thus, considering the 
averages ( iVFnn'l 2 ), the scattering cross section for 55° 
should present almost no shoulder around 4 eV, since in- 
terband excitations 41 and 42 are weak for that angle. 
Likewise, we expect the scattering cross section for 70° 
to have a weaker peak around 3 eV, because of the negli- 
gible contribution of excitations for nn' = 21, 22. We can 
appreciate these effects in Fig. [7[ Thus, given the peaks 
of the interband densities of Stoner states, the shift in 
the Stoner excitation maximum with varying scattering 
angle obeys to which interband excitations receive more 
weight, according to the corresponding matrix elements 
averages. This, in turn, depends on which umklapp pro- 
cesses gain more importance. Fig. [t] also shows that the 
maximum shift causes the broadness of the spectrum to 
increase with scattering angle. 

Turning to spin waves, it is interesting to note that the 
strength of the higher optical mode decreases with scat- 
tering angle. Thus it appears that umklapp scattering is 
unable to transmit to it sufficient oscillator strength at 
higher scattering angles. On the other hand, through 
comparison with the noninteracting cross sections we 
found that the lower lying spin wave seems less affected 
by scattering angle. 
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IV. DISCUSSION 

We wish to discuss some of the issues considered in 
this work pertaining to other theoretical and experimen- 
tal findings. First of all, as we have shown (cf. Fig. ^(b)), 
the maximum in the SPEELS spectrum does not neces- 
sarily correspond to the exchange splitting energy of the 
ferromagnet, since the peak is at 3 eV and A = 2 eV 
in our model. Thus, the common assumption that this 
peak allows us to estimate the magnetic moment [ p2"|]23"| 
through its correlation with exchange splitting should be 
reconsidered, both in bulk and surface studies. 

Also, the broadness of the spectrum is generally asso- 
ciated with a non constant exchange splitting over the 
Brillouin zone. While we agree a non constant exchange 
splitting will have this effect, we have seen that a most 
important source of broadness is umklapp scattering, to- 
gether with the structure of the interband densities of 
Stoner states of the material. In particular, even for 
specular or near specular scattering, will a model with 
rigidly split bands present a relatively broad spectrum. 
Still, according to our interband densities of Stoner states 
plot, for rigidly split bands we expect to see a maximum 
at exchange splitting. Unfortunately there are no reliable 
results for specular scattering in the case of Fe. p{| In 
this regard, the case of Ni would be very interesting to 
investigate in more detail. To begin with, the results of 
Kirschner, Rebenstorff, and Ibach and of Abraham 
and Hopster |l3[ ] are contradictory. Indeed, the former 
reported a broad maximum around exchange splitting 
for near specular scattering, while the latter stated that 
they see no sharp feature at that energy. According to 
our picture, Abraham and Hopster's result would be ex- 
plained only if the matrix elements W nn > are always weak 
for energies near exchange splitting. An accurate calcu- 
lation of the interband densities of Stoner states and of 
the weights W nn > would be most clarifying in this re- 
spect. For off specular scattering, however, both groups 
reported a weak dependence on scattering angle. This is 
plausible, according to our results. In Fig. ^|we see that 
the (lll^nn'l 2 ) may remain relatively unchanged for cer- 
tain scattering angle intervals, with, consequently, little 
change in the SPEELS spectra. 

As mentioned in the introduction, moreover, Abraham 
and Hopster interpret the onset of Stoner excitations 
found in their work in terms of the 3d band structure of 
Ni. Again, this could readily be verified having at hand 
the interband densities of Stoner states for this material. 
Indeed, these authors consider in particular interband ex- 
citations corresponding to nn' — 55, the onset of which, 
if they exist, could be easily identified in a graph like 
that in Fig. ||(a). A point still to be verified would be if 
the matrix elements for such excitations are sufficiently 
important. 

Spin waves, both acoustic and optical, have long been 
predicted in itinerant ferromagnets and subsequently ob- 
served through neutron scattering. Jl],£|,§| The question 



is, then, if these are observable in SPEELS. The question 
to be studied in the future is if there can be enough cou- 
pling between the incoming electron and those in the solid 
to excite an optical spin wave. In our model, it is umk- 
lapp processes that provide optical spin waves with the 
necessary oscillator strength to contribute significantly 
to the spectra. It could be, however, the acoustic waves, 
when present, drain most of the oscillator strength. This 
is plausible because it is known that acoustic spin waves 
in Fe arise upon hybridization of d electrons with sp elec- 
trons. Matrix elements with sp hybridization will be 
more important because of the larger s wave functions. 
This implies, of course, as our model does, that optical 
modes are mainly d in character. The case of Ni appears 
again to be different, since even a pure d band model of 
Ni shows acoustic spin waves. || 

The analysis presented in this work can prove useful 
more broadly in the understanding of ferromagnetism. 
Recently, Hirsch has presented a model of ferromag- 
netism without exchange splitting, in which spin polar- 
ization arises upon broadening of the spin-up bands rela- 
tive to spin-down. p0| ] If this mechanism plays an impor- 
tant role in itinerant ferromagnets, then the interband 
densities of Stoner states will change considerably be- 
cause other pairs of bands will be involved than those in 
the Stoner picture of ferromagnetism. Consequently, the 
predicted exchange scattering spectra will be different in 
both pictures. Thus SPEELS can prove a useful tool to 
validate or disprove Hirsch's model, possibly improving 
resolution previously. 

Finally, we would like to comment on the bulk vs. sur- 
face question. Our calculations here have focused on 
bulk properties. However, some authors have presented 
SPEELS as a technique more appropriate for surface 
studies. |2j|,|25| The reason is concern regarding the mean 
free path of electrons at the energies used in SPEELS. 
The question raised, however, is not simple and requires 
more detailed consideration, both theoretically and ex- 
perimentally. Still, most of the observations to date have 
been discussed in the light of bulk calculations. ]l2] , p^| , p0[ 
In this regard, a recent report on the electron dynamics at 
the surfaces of noble metals (Ag and Cu) is appropriate 
to mention. Burgi et al. J4l| have found that the dy- 
namics of hot electrons at surfaces can be dominated by 
bulk electrons. This offers more support to the premise 
that our results offer a sensible explanation for SPEELS 
results. 



V. SUMMARY AND CONCLUSIONS 

In this work we address the problem of the interpreta- 
tion of the SPEELS spectrum of itinerant ferromagnets. 
We find that considerably more information can be drawn 
from these measurements than has been recognized un- 
til now. We have found that the peaks of the spectra 
in the Stoner region are the image of a few interband 







densities of Stoner states of the material. These are very 
sensitive to the electronic structure of the material and il- 
lustrate very clearly the allowed and forbidden interband 
excitations. The important band couples in SPEELS are 
determined by the average weight of the squared matrix 
elements of the pendant Stoner excitations. In this re- 
spect, umklapp processes play a most fundamental role. 
Our model also predicts that optical spin waves should 
be excited in SPEELS experiments, with umklapp scat- 
tering providing the necessary oscillator strength. Our 
results allow us also to explain several of the features ob- 
served in SPEELS spectra. From the theoretical point 
of view, ab initio calculations of the interband densities 
of Stoner states and matrix elements W nn i would pro- 
vide closer look to the details of the mechanisms behind 
SPEELS. The differences between different ferromagnets, 
like Fe and Ni, could also be better understood. From the 
experimental point of view, measurements with higher 
resolution would be desirable, particularly for the detec- 
tion of spin waves. We think our results provide a good 
starting point for those further studies. 
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